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Abstract 

Let p n denote the probability that a random instance of the stable roommates 
problem of size n admits a solution. We derive an explicit formula for p n and compute 
exact values of p n for n < 12. 


1 Introduction 

Matching under preferences is a topic of great practical importance, deep mathemat¬ 
ical structure, and elegant algorithmics flU Q. A paradigmatic example is the stable 
roommates problem (3). Consider an even number n of participants. Each of the par¬ 
ticipants ranks all the others in strict order of preference. A matching is a set of n/2 
disjoint pairs of participants. A matching is stable if there is no pair of unmatched par¬ 
ticipants who both prefer each other to their partner in the matching. Such a pair is 
said to block the matching. The stable roommates problem is to find a stable matching. 
The name originates from the problem to assign students to the double bedroomes of a 
dormitory. Another application is the formation of cockpit crews from a pool of pilots. 

An instance of the stable roommates problem is defined by a preference table, in 
which each participant ranks all other n — 1 participants, most preferred first. For 
technical reasons we will assume that each participant puts himself at the very end of 
his preference list. Here are two examples for n = 4: 


1 : 

4 2 3 1 

1 : 

3 2 4 1 

(A > 3: 

3 4 12 

14 2 3 

(B) V: 

1 3 4 2 (1) 

2 14 3 

4 : 

3 2 14 

4 : 

12 3 4 

In (A), the marked matching (1,2)(3,4) 

is stable. In (B), 

there is no stable matching: 


whoever is matched with 4 can always form a blocking pair with someone else. Example 
(B) illustrates the fact that not all instances of the stable roommates problem have a 
solution. Let p n denote the probabilty that a random instance, chosen uniformely from 
all possible instances of size n, admits a solution. Our examples shows that 0 < p 4 < 1. 
The exact value is p 4 = 26/27. It has been computed by Pittel fU more than 20 years 
ago. No other values of p n are known exactly. Numerical simulations ||5l suggest that 
p n is a monotonically decreasing function of n that asymptotically decays like n -1 / 4 . 

In this paper we derive an explicit formula for p n that we use to compute exact 
values of p n for n < 12. And we discuss a generalization of this approach for odd values 
of n. 
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2 Stable Permutations 


A matching of size n can be interpreted as a permutation n of {1 ,n} that is com¬ 
pletely composed of 2-cyles. An obvious generalization is to allow arbitrary permuta¬ 
tions 7 r, but for that one needs to extend the definition of stability. A permutation 7r is 
called stable if it satisfies the two following conditions: 


V* : i does not prefer ir(i) to it 1 (i) 
i prefers j to 7r(i) => j prefers -rr(j) to i 


(2a) 


(2b) 


This definition includes permutations with fixed points. This is the reason why we've 
added each participant to the very end of his own preference list. But note that d2bl ) 
rules out that a stable permutation can have more than one fixed point. 

For permutations composed of 2-cycles (matchings) condition ( |2a| l is trivially satis¬ 
fied and condition ( |2b| l reduces to the usual "no blocking pairs" condition. Condition 
( f2al ) enforces each cycle of length > 3 to have a monotonic rank ordering: every mem¬ 
ber i prefers his predecessor 7r” 1 (i) to his successor 7r(i), and condition ( |2b| l prevents 
any member of the cycle to leave the cycle. 

The significance of stable permutations for the stable roommates problem arises 
from the following facts, proven by Tan (6j: 

1. Each instance of the stable roommates problem admits at least one stable permu¬ 


tation. 


2. If 7T is a stable permutation for a roommates instance that contains a cycle C = 
(vi,v 2 , ■ ■ ■ ,V 2 m) of even length, we can get two different stable permutations by 
replacing C by the 2-cycles (m, v 2 ), ■■■, (v 2m -i, v 2m ) or by (v 2 , v 3 ), ..., (n 2m , ig). 

3. If C is an odd-length cycle in one stable permutation for a given roommates in¬ 
stance, then C is a cycle in all stable permutations for that instance. 

These facts establish the cycle type of stable permutations as certificate for the exis¬ 
tence of a stable matching: An instance of the stable roommates problem is solvable if 
and only if the instance admits a stable permutation with no odd cycles. 

Consider again the two examples from the previous section. One can easily check 
that the permutation (1,2,3)(4) is a stable permutation for ( B ). Since it contains the 
odd cycle (1,2,3), (B) admits no stable matching. The permutation (1,3,4,2) is stable 
for (A). According to fact [2] its 4-cycle can be replaced by (1,3) (4,2) or by (3,4) (1,2), 
which are in fact both stable matchings. 

3 A Formula for p n 

The facts proven by Tan allow us to derive an explicit formula for the probability p n . 
The underlying ideas have already been discussed more or less in 0), but the formulas 
([8]) and ( | 1 3[ ) haven't been published before. We start with an integral representation for 
the probability that a permutation n is stable. 

Proposition 3.1. Let n be a permutation of {1,..., n} and let = {i : i = 7r(i)} denote 
the fixed points and M n = {i : i r(i) = 7r _1 (i) ^ i} the elements in two cycles of n. The 
probability that n is a stable permutation for a random instance of the stable roommates 
problem is given by 



( 3 ) 
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where integration is over the n-dimensional unit cube and 


Ar = {(i,j) i = 7r(j) v j = tt(*)} 


( 4 ) 


is the set of pairs of elements that are cyclic neighbors in n. 


Proof. A random instance of the stable roommates problem can be generated as follows: 
Introduce an n x (n — 1) array of independent random variable X i:j (1 <i^j< n), each 
uniformly distributed in [0, 1]. Each agent i ranks the agents j / % on his preference 
list in increasing order of the variables Xij. Obviously such an ordering is uniform for 
every i, and the orderings by different members are independent. The fact that each 
agent is at the very end of his hown preference list is taken into account by adding 
variables X rj = 1 to the set of random variables. 

Let P{tt\x, y) denote the conditional probability that the permutation 7r is stable given 
-X’in-(i) = Xi and X i7r -i W = y it and let F„ = {i : i = i r(i)} and M n = {i : tt(*) = 7r _1 (i) ^ i} 
denote the fixed points and two cycles of n. Then ( ]2aP tells us 


P(n\x,y) oc <d(xi~yi) 


ieM^UF,, 


(5) 


where 0 is the step function 


0 ( 2 ) 


1 z>0 
0 z<0 


and 8(z ) is the Dirac delta function. 

The second condition < [2b| ) is violated if X (? < and X 3i < x 3 for some (j,j) ^ D 7T . 
This happens with probability x t x 3 , hence 


P{x\x, y) oc ]^[ (1 -XjXi), 


( 6 ) 


which does not depend on y. 

Integrating ([5]) over y, gives a factor x, if i is an element of cycle of length three or 
more, a factor 1 otherwise. Adding the product Y\, eFr: 8(xi — 1) to ensure the constraints 
Xu = 1 finally allows us to integrate over the x,'s to obtain (3). □ 

Note that ([3]) differs slightly from the integral representation in [4j: Our integral 
is valid for any permutation n. If tt contains more than one fixed point, the integrand 
vanishes since the 5-function forces at least one of the factors in the product \\il-XiXj) 
to be zero and P(n) = 0 as it should. 

Obviously P( tt) depends on n only through the cycle type of n. Let a* denote the 
number of cycles of length k in n. We use the notation a = [l ai , 2“ 2 ,...] to denote the 
cycle type, including only those terms with au > 0. For n = 4, the only non-zero integrals 
are 

r 1 233 

A[ 2 2 ]) = J d A x{^-XiX 3 )(l-x 1 x i ){l-x 2 x ? ,){l-X 2 X i ) = — (7a) 

f 1 25 

-Pd 41 ]) = J d 4 z (1 - £1X3) (1 - X!X 2 x 3 X 4 = ( 7 b) 

r 1 1 

P([l 1 3 1 ]) = J d 3 x (1 — X\)(l — x 2 )(l — xs)x\X 2 X 3 = ■ (7c) 

Note that in the last integral, we have already done the trivial integration over S(x 4 — 1). 


3 



Proposition 3.2. Letp n (n even) be the probability that a random instance of the stable 
roommates problem has a solution. Then 

Pn= Y (-l) efa) c(a)P(a) > (8) 


where £ rl is the set of all cycle types of size n with even cycles only. The exponent e(a) 
is the number of even cycles of length > 4 in a, e(a) = J2k= 4 6 “/=• The factor c(a) is 
the number of permutations with cycle type a. 


c(a) 


n\ 


( 9 ) 


Proof. A matching of size n has cycle structure a = [2 n / 2 ], and there are (n— 1)!! match¬ 
ings of size n. Boole's inequality (aka union bound) then tells us that 


Pn < (n— 1)!! P([2"/ 2 ]), 


( 10 ) 


where equality holds if and only if the stability of different matchings were independent. 
This is not true in our case. Fact [2] from above tells us that stable matchings may come 
in pairs. Every stable permutation that consists of exactly one even length cycle of size 
z > 4 and (n — z)/ 2 cycles of size 2 corresponds to two stable matchings. These pairs 
have been counted twice in the sum in ( |1 0| ). The number of permutations of cycle type 
[2(™-z)/ 2 2 1 ] is n! ((n - z)\\ z) 1 and we get 

ra | 

Pn> {n— 1)!! P([2"/ 2 ]) — Y, ( _ V, P(l* n ~ z)/2 z1 ]) ■ (11) 

The > is again a consequence of Boole’s inequality. Equality in ( | 11P would only hold if 
the stability of pairs of permutations were independent events, but we know from fact[2] 
that stable pairs again may come in pairs: we have a quartet of stable permutation for 
each permutation that is composed of precisely two cycles of length > 4 and 2-cycles. 
Again we can express the corrections by P([]) and a combinatorial prefactor. Iterating 
this reasoning (which is of course the well known inclusion-exclusion principle) yields 

©• 

The formula © for the number of permutations of a given cycle type is well known. 
Yet we will give a short proof for completeness. Write down the cycle structure in terms 
of a fc pairs of parentheses enclosing k dots, like 


(•••)(•••)(••)(•) ( 12 ) 

for n = 9 and a = [l 1 , 2 1 , 3 2 ]. Now imagine that the n dots are replaced left to right with 

a permutation of {1,...Then the parentheses induces the desired cycle structure 
on this permutation. There are n\ permutations, but some of them result in the same 
"cycled" permutations. First, a cycle of length k can have k different leftmost values 
in (•••), which gives a factor k ak of overcounting. And pairs of parentheses that hold 
the same number of dots can be arranged in any order, which gives a factor afc! of 
overcounting. This yields © . □ 

Corollary 3.3. Let O n denote the set of all cycle types of size n that contain at most 
one fixed point and at least one odd cycle. Then 

l~Pn = Y (- 1 ) e(a) c( a ) J P( a ) , (13) 

a GO„ 
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Proof. Since P( a) = 0 if a has more than one fixed point, we can extend the sum to 
run over all cycle types with at least one odd cycle. Then the right hand side of fl!3| l is 
the probability that a random instance of the stable roommates problem has a stable 
permutation with at least one odd cycle. But this equals the probability that a random 
instance of the stable roommates problem has no solution. □ 


4 Evaluation of p n 


We already know the values of the integrals P(a) for n = 4, see 0. When we insert 
these values into 0 or fl3| we get 


p 4 = 3P([2 2 ])-6P([4 1 ]) 

l-p 4 = 8P([l 1 ,3 1 ])=4> 


26 

27 


0.962962... 


(14) 


the value computed by Pittel in 1993 I®. It seems straightforward to compute p n for 
larger values of n, since all we need to do is to evaluate and sum the corresponding 
integrals P(a). This is not easy, however. Pittel wrote "For n = 6, the computations 
by hand become considerably lengthier and we gave up after a couple of half-hearted 
attempts." 

The computations become "lengthier" for two reasons: the number of integrals in 
0 and ( |13| ) increase with n, and the evaluation of each individual integral gets harder. 
Let us first look at the number of integrals: 


Lemma 4.1. Let p{n) denote the number of unordered partitions of n, and let n be 
even. Then 

|£n| = p(|) (15a) 

\O n \ =p(n)-p(n-2)-p^j . (15b) 

Proof. From ka k = n, or from glancing at m it is obvious that there is a one-to-one 
correspondence between the set of cycle types of size n and the set of integer partitions 
of n. 

Every cycle type a s £„ corresponds to a partition of n into even numbers and vice 
versa. Every partition of n into even numbers corresponds to a unique partition of n/2 
and vice versa—simply divide or mutiply all parts of the partition by two. This proves 
( fl5al >. 

The number of all cycle types is p(n), and the number of all cycle types that contain 
at least two fixed points is p(n — 2). Hence the number of cycle types that contain at 
most one fixed point is p{n) — p(n — 2). For \O n \ we also need to subtract the number of 
cycle types with even cycles only, which is p(n/ 2). This proves d 1 5b[ >. □ 


There is no closed formula for the partition numbers pin), but they are known for all 
n < 10 000 0. And we need p{n) only for small values of n to get 


n 

4 

6 

8 

10 

12 

14 

16 

18 

\£n\ 

2 

3 

5 

7 

11 

15 

22 

30 

\Ou\ 

1 

3 

6 

13 

24 

43 

74 

124 


In this regime of n, the number of integrals is no problem. So let us turn our attention 
to the individual integrals. 

When we expand the product in 0, we get a sum of easy-to-integrate terms of the 
form x J 1 • x’p-, but there too many terms to be integrated by hand. 
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Pi 

©: 0.20 sec. 

( fl3l l: 0.02 sec. 


P6 P8 P 10 

19.8 sec. 5 min. 20 min. 
3.5 sec. 6 min. 25 min. 


P 12 

15.5 days 
13.9 days 


Table 1: Times to compute p n according © or 113 - 


Lemma 4.2. A full expanion of the integrand in © yields terms , where 

/(a) = -n{n — 3) + ag + a 2 . (16) 

Proof. If we expand the integrand, each factor in the product 


n (i -x^) d7) 

i<3 

doubles the number of terms. Hence we need to show that ( | 1 6[ ) is the number of factors 
in this product. Think of the n variables ag as the vertices of a graph G. Each factor 
(1 — XiXj) in dz) corresponds to an edge of G. Without the constraint (i,j) ^ D n , G is 
the complete graph with |n(n — 1) edges. Each cycle of length k > 3 in a corresponds 
to a cycle in G with k edges that are removed from the complete graph. Each cycle of 
length 2 corresponds to an edge that is also removed. This gets us 

/(a) = ^n(n - 1) - ^ ka k - a 2 = \ n ( n - 1) - (X! kak ~ 2 “ 2 - ai ) - ° 2 
k >3 V k ) 


and fl6l l follows from J2k ka k = n - 


□ 


The maximum number of terms arises for pure matchings, i.e., for a 2 = n/2 and 
a\ = 0. It reads 2 4 , 2 12 , 2 24 , 2 40 and 2 60 for n = 4,6,8,10 and 12. Hence it is no surprise 
that Pittel gave up on the integrals for n = 6. The integration is better left to a computer. 

We used the computer-algebra system Mathematica ® for the exact evaluation of 
the integrals P(a). Figure [l] shows the Mathematica code that sets up the integrand 
and performs the integration. The full Mathematica code is available online |9). 

Using our Mathematica code, we computed the values ofp n for n < 12 both from © 
and (as a crosscheck) from l |13| l. The results are 


Pe 


Ps 

P 10 


181431847 
194400000 
809419574956627 
889426440000000 
25365465754520943457921774207 
28460490127321448448000000000 


= 0.93329139403292181070... 


= 0.91004667564933981499... 


= 0.89125189485653484085... 


13544124829485098788469430650439043569062157071 
“ 15469783933925839494793980316271247360000000000 
= 0.87552126696367780620.... 


(18a) 

(18b) 

(18c) 

(18d) 


The values of the corresponding individual integrals are listed in Tables [2] an d0 

We ran our Mathematica code on a computer eguipped with 2 Intel® Xeon® CPUs 
E5-1620 with 3.60 GHz clock rate and 32 GByte of memory. The total computation 
times are shown in Table [l] Table [3] also shows the times to compute the individual 
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Integrand[a ] := Module!(* computes integrand corresponding to cycle pattern a *) 
{n,i,j,l,result,cycle}, 

If[a[[1]]>1,result=0, (* more than one fixed point *) 
n=Sum[k*a[[k]],{k,1,Length[a]}]; 

If[a[[1]]>0, (* take care of fixed point *) 
n=n-1;result=Product[{l-x[i]),{i,l,n}], 
result = 1 (* no fixed point *) 

]; 

result=result*Product[(l-x[i]*x[j]),{i,l,n-l},{j,i+1,n}]; 

(* remove 2-cycles from product *) 

result=result/Product[(l-x[2*i-1]*x[2*i]),{i,l,a[[2]]}]; 

(* cycles larger than 2 *) 

result=result*Product[x[i],{i,2*a[[2]]+l,n}]; 

For[cycle=3,cycle<=Length[a],cycle++, 
l=Sum[i*a[[i]],{i,2,cycle-1}]+1; 

For[i=l,i<=l+cycle*(a[[cycle]]-1),i+=cycle, 

For[j =0,j<cycle,j++,result=result/!l-x[i+j]*x[i+Mod[(j+1),cycle]])] 

] 

] 

]; 

result 


P[a ] := Module] 

{y,n,k}, 

n=Sum[k*a[[k]],{k,1,Length[a]}]; 

If[a[[1]]>0,n=n-1]; 
y = Integrand[a]; 

For[k=n,k>=l,k--,y=Integrate[y,{x[k],0,1}]]; 

y 

]; 


Figure 1: Mathematica code to compute the integrals P( a) The procedure 

Integrand [a] returns the integrand as a function of variables x [ 1],... ,x [n] (or x[n-1] 
if the cycle type a contains a fixed point), the procedure P [ a ] evaluates the integral by 
exactly integrating variable by variable. 
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a 

P ( a ) 

a 

P(a) 

[2 3 ] 

448035973 


38077 

5832000000 

86400000 

[2M 1 ] 

307841 

[O.S 1 ] 

26257 

144000000 

777600000 

[fill 

2591729 

[921 

1742111 

1° J 

11664000000 

l' 3 J 

7776000000 

19 4 1 

1245959394495647 

[l 1 ,? 1 ] 

49958102093 


107585022182400000 

384232222080000000 

[2 2 ,4 1 ] 

5211637894488503 

[O^s 1 ] 

441974732789 

26896255545600000000 

12807740736000000 

[2S6 1 ] 

914248620325799 

[iSaM 1 ] 

1249592153 

53792511091200000000 

9605805552000000 

14 2 1 

1493807915753 


58105985423 

l 4 J 

1195389135360000000 

25615481472000000 

[a 1 ] 

622186155317 

[2\3 2 ] 

76670733315619 

1° J 

498078806400000000 

4482709257600000000 



P 1 ^ 1 ] 

58105985423 



25615481472000000 

[2 5 ] 

433857166916418660757431885203 


1882697003227025150390719 

322741958043825225400320000000000 

819662115666857715302400000000 

[2 3 ,4 1 ] 

4794693488032751578104859937 

[i 1 ,3 3 ] 

158398327239405983477 

322741958043825225400320000000000 

512288822291786072064000000000 

[2 2 ,6!] 

726158117631681830112186713 


2765878679393466620633 

645483916087650450800640000000000 

409831057833428857651200000000 

ol A 2 

1 94089601969271248978571831 

[iSz 2 ^ 1 ] 

4336602947669955694769 

Z , 4 _ 

1 1290967832175300901601280000000000 

32786484626674308612096000000 

[2S8 1 ] 

18812621042800384360939621 


126601947989502609349 

258193566435060180320256000000000 

409831057833428857651200000000 

[4 1 ,6 1 ] 

10678226865621944175135083 


633196266619396193087 

2581935664350601803202560000000000 

2049155289167144288256000000000 

[10 1 ] 

42708804188035567140443357 


789921304062168675601 

10327742657402407212810240000000000 

117094587952408245043200000000 



[O.S 1 ] 

1265995491264426770353 



4098310578334288576512000000000 



[2 2 ,3 2 ] 

2918990176269285877130918549 



2581935664350601803202560000000000 




18845369089082632479619357 



258193566435060180320256000000000 



[3 2 ,4 4 ] 

21410287713579117222366871 



5163871328701203606405120000000000 



[3 1 ,7 1 ] 

610894828667022260751797 



147539180820034388754432000000000 



^21 

8541874436295301342281403 


2065548531480481442562048000000000 


Table 2: Probabilities P(a) for n = 6,8,10 (top to bottom). Cycle types with (right) and 
without (left) odd cycles. 


integrals for n = 12. Some of theses integrals (marked with a *) could not be computed 
by the simple iterative scheme in Figure [T] because Mathematica ran out of memory. 
In these cases we expanded the integrand in a polynomial in the variable x n (or x n -i 
if there is a fixed point) and applied interative integration to each coefficient of this 
polynomial. This reduces the memory consumption, but it slows down the computation. 
With a larger memory (like 64 GByte instead of 32 GByte), this could have been avoided 
and pi 2 could have been computed somewhat faster. 


5 Odd values of n 

For odd values of n there are no stable matchings, of course. But there are still sta¬ 
ble permutations: Tan’s results listed in Section [2] as well as Proposition |3. 1 1 also hold 
for odd values of n. This allows us to generalize the stable roommates problem to odd 
values of n. The most obvious generalization is to accept one fixed point, i.e., to reject 
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a 

P(a) 

Time [sec.] 

[ 2 6 ] 

325899908494883644126440199857602193757211429627 

2572934463890545624774134806202233860915200000000000 

265018 

[ 2 M 1 ] 

1209115974791734652605681563324122140963407221 

1225206887566926487987683241048682790912000000000000 

265 091 

[ 2 2 ,4 2 ] 

16288072152327610053000950409164225186650151 

4288224106484242707956891343670389768192000000000000 

205 089 

[4 3 ] 

231703173390597300042186053017177445722753 

25729344638905456247741348062022338609152000000000000 

206115 

P 3 ^ 1 ] 

3378294177941932509053172872298924486923764591 

51458689277810912495482696124044677218304000000000000 

49 493 

PP^e 1 ] 

417880592074077264470531487240272595070729 

2144112053242121353978445671835194884096000000000000 

49 220 

[ 6 2 ] 

1853330912748299530044034784734880537462353 

205834757111243649981930784496178708873216000000000000 

49 293 

pps 1 ] 

508893194633666952579907861671385135829521 

134007003327632584623652854489699680256000000000000 

47 303 

[4S8 1 ] 

4412925241742005785167715449536219676971 

490082755026770595195073296419473116364800000000000 

47 520 

[ 2 1 , 10 1 ] 

53484730261191253361608747405394814514357 

274446342814991533309241045994904945164288000000000 

53 836 

[ 12 1 ] 

14826641669894164340076941832557808383893 

1646678056889949199855446275969429670985728000000000 

53299 

[ 1 P 2 1 .3 1 ] 

122503966894472107602242737308438169403 

920820472257490446118689304539955200000000000 

92 605 

[1*,2 1 ,3 3 ] 

1222543880169122622560877738575059 

93543667022983156431104945223106560000000000 

42 996 

- 1 ) 2 2 j 3 1 ,4 1 ] 

193959334006722457965074605079586629657 

618791357357033579791759212650849894400000000000 

7 038 

[l\3\4 2 ] 

2855025200767172513735081732106863 

5729549605157718331405177894915276800000000000 

8 354 

[ 1 \ 2 3 ,5 1 ] 

8437895290055710585350317910566101247317 

1237582714714067159583518425301699788800000000000 

10 784 


176230039945164631684723423501203937 

353595061346876331309576692943342796800000000000 

10498 

1 5 2 1 ,4 1 , 5 1 ] 

1795748861495201189078302039999290739 

137509190523785239953724269477966643200000000000 

9 092 

1 5 2 1 ,3 1 , 6 1 ] 

2694291584800385097759305707384840499 

206263785785677859930586404216949964800000000000 

8314 


2174740904996317920876887889334183 

4365371127739213966784897443744972800000000000 

8 376 


1695852720842076492466118915028628181 

5412169306329739764942500402194022400000000000 

8 059 


88077935438211707375963113857429259 

176797530673438165654788346471671398400000000000 

8 082 


19270992061340957619880212582236803 

38674459834814598736984950790678118400000000000 

8 098 


1197147888812853403164542654655246797 

91672793682523493302482846318644428800000000000 

8 096 

[iMi 1 ] 

3699232196202777873480674898554033731 

7425496288284402957501110551810198732800000000000 

8001 

[ 2 3 ,3 2 ] 

161499154693883709213457621140881273357670713 

2450413775133852975975366482097365581824000000000000 

209 450 

p 2 ,3\5p 

4348643545825433788892694617856351275828603 

1143526428395798055455171024978770604851200000000000 

56207 

PPaM 1 ] 

357327697339220191285941924523831564051 

1829642285433276888728273639966032967761920000000 

212 763 

[3 4 ] 

9836824308843120655187019024812769030613 

1089072788948379100433496214265495814144000000000000 

216063 

\2\Z\1 1 ] 

91054335946285045516721350625722458874631 

466745480977876757328641234685212491776000000000000 

50230 

[S 2 ^ 1 ] 

618748213165565813756023362735829302961763 

68611585703747883327310261498726236291072000000000000 

49 614 

[ 2 1 ,5 2 ] 

26742627021755978677974945561880197844453 

137223171407495766654620522997452472582144000000000 

59 920 

PP^s 1 ] 

61829720130534204153121031185033846815523 

6861158570374788332731026149872623629107200000000000 

59 587 

[3 1 ,9 1 ] 

20608745217756332777239346633128265563123 

2287052856791596110910342049957541209702400000000000 

53 642 

[5 1 ,? 1 ] 

17650822382212529975478949518938933533441 

1960331020107082380780293185677892465459200000000000 

53495 


Table 3: Probabilities P( a) for n = 12 and the times to compute them. Times marked 
with * refer to a slower, more memory efficient integration procedure (see text). 
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one participant from the dormitory (or put him into a single bedroom), and to ask for a 
stable matching of the remaining n — 1 participants. Let p n (for n odd) denote the prob¬ 
ability that a random instance admits such a solution. Following the same reasoning as 
in Proposotion |3.2| and Corollary |3.3| we get 

Pn = E (-!) e(a) c( a )P( a ) , (19) 

aefi 

1-Pn= 53 (-!) e(a)c ( a ) p ( a ) > (2°) 

a £03 


where £ 3 is the set of all cycle types of size n consisting of one fixed point and even 
cycles and O 3 is the set of all cycle types of size n that contain at least one cycle of odd 
length > 3. Table [4] lists the values of the corresponding integrals P( a) for odd n < 11. 
The resulting values of p n are 


P 3 


^ = 0.75 
4 


P 5 

P7 

P9 


Pll 


4075 


= 0.5895543981481481... 


6912 
246462083 


= 0.4754284008487654... 


518400000 

11365049284140796201 
29144725585920000000 
176967745750762518431538515329 
546441410444571810201600000000 


= 0.38995218021992023... 


= 0.3238549318705289 ... 


(21a) 

(21b) 

(21c) 

(21d) 

(21e) 


It seems counterintuitive that p 2 k -i < but note that the enforced fixed-point for an 
odd number of participants represents someone who is happy to be matched with any¬ 
body else. This high destabilizing potential is a result of the rule that every participant 
has to put himself at the very end of his preference list. 


6 Conclusions and Outlook 

We have seen that p n , the probabilty of a random instance of the stable roommmates 
problem of size n to admit a solution, can be expressed as a sum over cycle types of 
permutations of size n. Each term in the sum is an integral with an exponential number 
of terms. The latter restricts an exact evaluation of p n to n < 12. In spite of this 
limitation, the method is far more efficient than the exhaustive enumeration over the 
[(n — l)!] n_1 different instances of size n. For n = 12, this number is 4.1 x 10 83 , or 4100 
times the number of atoms in the visible universe (which is usually estimated as 10 8 °). 

Our results for n < 12 don't shed new light on the ultimate behavior of p n as n 
becomes large, but they suggest that exact evaluation of p n for any larger values of n is 
likely to be infeasible without some unexpected new approach. 

The approach outlined in this paper can easily be modified to work for the stable 
matching problem on general graphs, where each participant corresponds to a vertex 
of a graph G and ranks only those participants adjacent to him in G. If G is the complete 
graph, we recover the stable roommates problem. In the case of bipartite graphs G 
(known as stable marriage problem) we have p n = 1. For non-bipartite graphs, p n 
seems to be a monotonically decreasing function of n that may or may not approach a 
non-zero value, depending on the number of short cycles in G m. 
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a 

P(a) 

a 

[1S2 1 ! 

1 

4 

[3 1 ] 

[1\2 2 ] 

833 

20736 

[^s 1 ] 

[l 1 ) 4 1 ] 

1 

2304 

[5 1 ] 

[1\2 3 ] 

110831617 

23328000000 

[2 2 ,3 1 ] 

[1S2S4 1 ] 

797731 

23328000000 

[^s 1 ] 

[1S6 1 ] 

6541 

2916000000 

[3M 1 ] 

[l\3 2 ] 

2183 

972000000 

[7 1 ] 

[l\2 4 ] 

12242957448855683129 

27541765678694400000000 

P 3 ^ 1 ] 

[l\2 2 ,4 4 ] 

2998148628185909 

1311512651366400000000 

[2 2 ,5 1 ] 

[1 1 , 4 2 ] 

184134811312313 

27541765678694400000000 



3617070987119831 

27541765678694400000000 

[3 3 ] 

[1S8 1 ] 

26303739761759 

3934537954099200000000 

[2 1 ,7 1 ] 

[l\2\3 2 ] 

1206877128048157 

9180588559564800000000 

p 1 ^ 1 ] 

[1S3S5 1 ] 

2455964944171 

367223542382592000000 

[4 1 ,5 1 ] 
[9 1 ] 

[l\2 5 ] 

88853486478784120344992170351 

2581935664350601803202560000000000 

I2 4 ,3 1 ] 

[l 1 , 2 3 , 4 1 ] 

4572509990406797552502341 

34425808858008024042700800000000 

[2 3 ,5 4 ] 

[l 1 , 2 1 , 4 2 ] 

268054718171660435931803 

860645221450200601067520000000000 



1168831786137020235667067 

172129044290040120213504000000000 

[2\3 3 ] 

[l 1 A 1 ,#] 

743715115011041403407 

57376348096680040071168000000000 

[2 2 , 7 1 ] 


100516822545753167453891 

322741958043825225400320000000000 


[l 1 ,10 1 ] 

66933419040890225282203 

5163871328701203606405120000000000 


[1\2 2 ,3 2 ] 

4386643900008678909343237 

645483916087650450800640000000000 

[3 2 > 5 1 ] 


531499853646597948383 

40983105783342885765120000000000 

[3 1 ,4 2 ] 

1 ,2 1 ,3 1 ,5 1 ] 

804392338445761377188767 

2581935664350601803202560000000000 

p 1 ,^] 

[1 4 ,5 2 ] 

16733378688533122315949 

1290967832175300901601280000000000 

[3S8 1 ] 

[l 1 ,3 1 ,7 1 ] 

19128779897378689455131 

1475391808200343887544320000000000 

[4 1 ,7 1 ] 

[ll 1 ] 

62675660640300931114214381 


309832279722072216384307200000000000 


P(a) 

1 

8 

491 

27648 

191 

82944 

5103637 

2592000000 

1945639 

9331200000 

336349 

18662400000 

558779 

31104000000 

39406434169244998649 

220334125429555200000000 

234360972607515209 

14688941695303680000000 

3502136387768779 

2937788339060736000000 

374799675933251 

4896313898434560000000 

12476274579169301 

10492101210931200000000 

16811008475015879 

220334125429555200000000 

671436255551711 

8813365017182208000000 

1864835590786319 

24481569492172800000000 

710107424563570828306588840739 

51638713287012036064051200000000000 

109197089334060411167876570117 

103277426574024072128102400000000000 

14352373021321999225705658471 

206554853148048144256204800000000000 

183002200715285406357445301 

45901078477344032056934400000000000 

204627732127480795488157591 

2950783616400687775088640000000000 

547557978971950371021494551 

137703235432032096170803200000000000 

2461017693717356460362113427 

619664559444144432768614400000000000 

41866526759821300816071211 

206554853148048144256204800000000000 

373490614662460067378083 

1844239760250429859430400000000000 

182277891278802756936253723 

45901078477344032056934400000000000 

62737571161936687651226813 

309832279722072216384307200000000000 

17908596195396917111551979 

88523508492020633252659200000000000 

6963996535691809265274221 

34425808858008024042700800000000000 


Table 4: Probabilities P( a) for n = 3,5,7,9,11 (top to bottom). 
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